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Capillary rise is one of the most well-known and vivid illustrations of capillarity, however, there is 


no complete solution has been obtained except some segmental solutions for asymptotic regimes. In 


this paper, we use singularity-free Bush’s equation and successfully obtain its series and perturbation 


solutions. The solutions reveal that the capillary rise dynamics is mainly controlled by the Bond 


number and the Galileo number, of which the Bond number is a key parameter on the existence of 


the solution. Due to the poor rate of convergence of the series solution, we propose an approximate 


pga 


analytic solution h(t) = 272252 fı h 


=i /i29 t 
(y 2) t4 at’ e ( te ) l where a is a parameter. 


The proposed analytic solution has been verified numerically. 
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INTRODUCTION 


Capillary rise is one of the most well-known and vivid 
illustrations of capillarity. Knowledge of capillarity laws 
is important in oil recovery, civil engineering, dyeing of 
textile fabrics, ink printing, and a variety of other fields. 
It is the capillarity that brings water to the upper layer of 
soils, drives sap in plants, or lays the basis for operation 
of pens [1-25, 27-33, 35-38]. The capillary action has 
been studied by numbers of masters in science. Leonar- 
do da Vinci (1452-1519) recorded the effect in his notes 
and proposed that mountain streams may result from 
capillary rise through a fine network of cracks. Jacques 
Rohault (1620-1675): erroneously suggested that capil- 
lary rise is due to suppression of air circulation in nar- 
row tube and creation of a vacuum. Geovanni Borelli 
(1608-1675) demonstrated experimentally that h ~ 1/a. 
Geminiano Montanari (1633-87) attributed circulation in 
plants to capillary rise. Francis Hauksbee (1700s) con- 
ducted an extensive series of capillary rise experiments re- 
ported by Newton in his Opticks but was left unattribut- 
ed. James Jurin (1684-1750) independently confirmed 
h ~ 1/a; hence Jurin’s Law. Albert Einstein, whose first 
paper was submitted to Annalen der Physik in 1900, was 
on capillarity [1]. 

Washburn [3] developed the equation to describe the 
rate of penetration of liquids into small cylindrical cap- 
illaries based the laws of hydraulics. He assumed that 
the distance traveling in the initial turbulence period is 
negligible, and the Poiseuille region covers practically the 


entire flow. According to Poiseuille’s law, neglecting the 
air resistance, the rate of flow through a cylindrical takes 
the following form: Q = TAE (at + Lsa?), where a is the 
radius of the tube, h is the capillary penetration height, 


L, is the slip length. Q is the rate of flow and for a 
cylindrical tube it is given by Q = ra?44. The total 
effective pressure AP, consists of the external pressure 
Pe, the hydrostatic pressure P, and the capillary pres- 
sure given by Young-Laplace equation AP = 222288 
Combining the above relations and neglecting the s- 
lip length gives the so-called Lucas-Washburn equation: 
dh (Pa + Pp + 22298°). Tf there is no external pres- 
sure, a commonly used form of Lucas-Washburn equation 
is given by 8uh 44 +pgha?siny = 2accos6, where y is the 
incline angle of the tube with respect to the horizontal 
surface. 
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FIG. 1: The incline capillary rise 


When the gravitational force is negligible, a well- 
known Washburn’s law can be obtained with initial con- 
dition h(0) = 0: h = \/*%52°*¢, which predicts burst-like 
behavior with the velocity being infinite at zero time. 


This singularity of the solution highlights a deep inconsis- 
tency of the above equation [20]. Rideal [4] also obtained 
this equation through a rather simple derivation. 


Taking into account of the momentum of the liquid in 
the tube and end-effect drag on the fluid entering the 
tube, Brittin [6] derived a more rigorous formulations of 
mee Washburn equation as follows: h4% dh + Š(%)? + 

zaz uh 
viscosity, ø is the liquid- air surface tension 0 = yzyv, 0 is 


+gh= Zo cos, where p is density, u is the 


the wetting angle of the liquid, h is the height of capillary 
rise, a is the capillary radius, and g is the acceleration of 
gravity. This equation assumes the Poiseuille flow profile 
throughout the capillary. 


The most papular equation of capillary rise dynamics 


is hae + (By + ooo uh Fe +gh= Z0 608 0. Although 


this equation has been widely used in on literature, how- 


ever it has a singularity at t = 0, namely h = un — oo 
as t + 0, which would lead to an ill-posed problem. This 
equation can not even deal with the natural initial con- 
ditions h(0) = h(0) 
h(0) = \/20 cos 6/(pa) [5], neglecting such a logical draw- 


= 0. A formal remedy is to take 


back as the acceleration of liquid front at zero time being 
infinite, more criticism on this equation can be found in 
(8, 20, 31]. 

This singularity problem was firstly pointed out by 
Szekely et al.[8], who successfully removed the singulari- 
ty problem by composing the correct energy balance for 

2 
the entry flow: (h+ fa) GP +1.225(%)?4 Si Eh% +gh= 
r cos 6. Xiao et al.[28] considered the entrance pressure 


loss effects by using Dreyer’s model [13] and dynamic 
contact angle effect by Newman’s model [7] to further 
modified a Washburn equation: n +1. 028a) $} + 
0.958( 42)? + 4 (h +0.25a) È + gh = 22 cos 6. Bush [36] 
derived a similar equation with a diferenti in the coeff- 


cient of (4)?. 
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For centuries of investigation on the capillary rise dy- 
namics, it is incredible that no complete solution has been 
obtained yet. Ones have to find some information for 
asymptotic regimes. How to match the different solu- 
tions between the regimes has never been solved and still 
remains as an unsolved problem. It would be natural 
thinking that ones can not satisfy those asymptotic solu- 
tions, and must be looking for a complete solution that 
can cover whole domain of time. In the following, we will 
mainly focus on Bush’s model and give its series solution 


and approximate analytic solution. 


To make the paper self-contained, the paper is or- 
ganised as follows. After this introduction, in Section 
2 the Bush’s equation has been formulated. In Section 
3 the perturbation solution and the series solution are 
obtained. In Section 4 an approximate analytic solution 
is proposed by using the initial conditions and the se- 
ries solution. In Section 5 the numerical verification and 
discussions have been presented. Finally, section 6 con- 


cludes the paper. 


BUSH’S FORMULATIONS OF CAPILLARY RISE 
DYNAMICS 


In physics, the capillary dynamics process is a strug- 
gle between surface tension with the combining effect of 
wall resistance and gravity. In the initial stage, the sur- 
face tension plays a leading role, and later the combining 
effect of wall resistance and gravity takes over. 

Capillary rise dynamics is governing by a combina- 
tion of gravity, viscosity, fluid inertia and dynamics pres- 
sure. From continuum mechanics, the dynamics pro- 
cess must satisfy the conservation law of momentum: 
tim (t) 22] = F + fg pv -v ndS, where the second 
term on the right-hand side is the total momentum flux 
ra?ph? = mh, so the force balance on the column capil- 
lary may be expressed as 


(pna°h+ ma )h= 2naccos§ —pra*hg 


aia ca R 
inertia added mass capillary force weight 
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viscous force 


(1) 


dynamical pressure 


We can estimate the added mass m, from the change in 
kinetic energy as column rise from h to h + Ah, namely 
AE; = 4A(mU?), where m = me+mo+Meo (mass in the 
column, in the spherical cap and all other mass, respec- 
tively). In the column, me = ma?hp, u = U; in the spher- 
ical cap, Mmo = 2z a3 p; and in the outer region, radical 
flow extends to oo, nut u(r) decays. Volume preservation 
U/2. Con- 


tinuity gives 27a?u,(a) = 2rr?°ur(r) and leads to u,-(r) = 


requires: ta?U = 2ra?ur (a) leads to u,(a) = 


az u,(a) = £U . Thus the kinetic | in the far field 
is MU? = $f ur(r)?d $f. ulr 
3 


which leads to mao = Vora. . Hence the change in ki- 


)2 p2rrdr, 
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netic energy will be 


1 1 
OE, = zA Me + Mo + Mmo )U? + zm2UAU 


1 1 
= 5 AmU" + 5 (Me + mo + Meo)2UAU 


1 
= 5 Ala pAh)U* + (ma*hp + T a8 p)UAU. (2) 


The added mass is estimated as Mg = Mmo + Mo = Tra’? p- 
For a typical capillary rise, a = 1074 ~ 5 x 1074m, u = 
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FIG. 2: The dynamics of capillary rise 


1.5 x 1073#8, p = 750kgm™? and top velocity can be 
at V = 0.04 ~ 0.778m/s, hence the typical capillary 
rise Reynolds number Re = eve = 10 ~ 195. At such 
low Reynolds number, if the tube is smooth one the flow 
in the tube is fully developed Poiseuille (laminar) flow, 
hence the viscous force is T, = —*#h, 

Combining all above relations, we have Bush equation 
of capillary rise dynamics [36]: 


7 @h 1 dh 


6” ap + ace 


8u , dh 
pa? dt 


(h+ 4 +gh= a (3) 
pa 

and initial height and velocity boundary conditions 

h(0) = 0, h(0) = 0. 

Although no complete solution of Eq.(3) has never 
been obtained, some segmental solutions (as shown in 
Figure 3) for asymptotic regimes have been proposed by 
reduced form of Eq.(3). (1) Initial regime: h ~ 0, h ~ 0, 
infers h(t) = 2 coos? hence h(t) = $ 759324. (2) Quasi- 
steady stage (Lucas-Washburn equation), where the cap- 


illary force is compensated by gravity and viscous drag, 
the asymptotic solution corresponding to the short-time 


limit (t — 0) is given by A(t) = ,/750°°t. (3) Low- 
viscosity limit (Quéré equation), the capillary rise is giv- 
20 cos 0 

pa `’ 
to stationary level, established by the balance of gravity 


en by h = (4) The capillary liquid will rise 
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FIG. 3: Asymptotic regimes of capillary rise dynamics [36] 


However, those segmental solutions have not been 
matched into a single solution that valid for whole do- 


main of time. 


SOLUTION OF BUSH’S EQUATION OF 
CAPILLARY RISE DYNAMICS 


Perturbation solution 


Bush’s equation (3) can also be transferred into a di- 
mensionless form as follows: 


dz 7 Bo dz 7 dz 


1l6cos@ dz 
74 tz—-1=0 
oat Beel ode 


T zZ T 
Bo-Ga_ dr 
(4) 


and the initial height and velocity boundary conditions 
are changed to z(0) = 0, (0) = 0. In which, z = 4, 
t, H = 229088, T = E, Bond number Bo = 4 
represents the ratio of gravity and surface tension, and 


: — Pga? 
Galileo number Ga = n2 


T = 


represents the ratio of gravity 
and viscosity. 


From the Eq.(4), we can notice that the Bond num- 
ber Bo is playing a vital role in capillary rise dynamics, 
because it stands for the ratio of gravity and surface ten- 
sion. In the very beginning (within the diffusion time 
2y surface tension is mainly compete with gravity, on- 
ly after the diffusion time the viscous drag will come to 
against the surface tension. In other words, from the be- 
ginning of capillary rise, the Bond number Bo enters the 
process and later on the Galileo number Ga joins. 


We choose the Bond number Bo as a perturbation 
parameter, the perturbation solution can take the for- 
m z(T) = Xp; Bo™” fa(T), which leads to a solution as 


follows 
1 6 1 ,42 
2(7T) = Bo!7 cos 0)r? — 3 (343 cos? @)r* 
1 144 34 
PLETE cos? 0)r 
1 3456 1 


= Bo! 1715 Ga cos* 0)r° + O(Bo~”) (5) 


Exact Taylor’s Series solution 


Eq.(3) is a nonlinear ordinary differential equation, up 
to now there is no closed form solution has been found- 
ed. Fortunately, it has an exact Taylor’s series solution 
h(t) = Xpo Cnt” as follows 


h(t) = 6a cos > 30 cos0(7Tgpa? + 240 cos 0) 
T pa? 343p7a° 
864107 cos? 0 
~ —-1715p3a7 


tt 
t + O(t®) (6) 


The correctness of this solution has been verified by using 
symbolic mathematics software Maple. Surprisely to see 
that the series solution is exact same as the perturbation 
solution in Eq.(5). The solution in Eq.(6) gives a im- 
portant result, namely, the non-zero initial acceleration 
h(0) = +2 2298 , which later will be used in the proposal 
of approximate analytic solution. With the above solu- 


tion in Eq.(6), there is no longer needs to divide the cap- 
illary dynamics process into several asymptotic regimes, 
since the solution in Eq.(6) is valid for the whore range of 
time domain. From Eqs.(5,6) we can see that the Bond 
number Bo is really vital for the existence of the solu- 
tion. For those equations without the added mass term 
Ta means that they have zero Bond number, which will 
lead to the infinite of height and velocity ! 

The academic value of the series solution is enable us 
to see the importance of the added mass that can void 
the singularity. However, the rate of convergence of this 
solution is very slow and useless for practical application. 
It would be natural attempts to propose a simple solution 
for the capillary rise. 


APPROXIMATE ANALYTICAL SOLUTION 


General speaking, the whole dynamics process can be 
qualitatively described as follows: in an infinite reser- 
voir, the capillary rise and velocity are zero in the initial 


state. Due to the effect of the surface tension, the cap- 
illary liquid obtains the initial acceleration (the initial 
acceleration must never be zero), begins to rise at a rel- 
atively uniform velocity, and the surface tension plays a 
dominant role in the ascending phase; however, as the 
capillary rises wall frictions and gravity begin to work 
in an attempt to prevent the rise of the capillaries, and 
their joint action succeeds in decelerating the capillaries 
at a point in time until the capillaries are finally stopped. 
Surface tension and wall resistance, gravity to achieve u- 
nity of opposites, capillary dynamics process are over, all 
attributed to calm. 

From the physical understanding of capillary rising 
phenomena, the capillary liquid will rise to stationary 
level, established by the balance of gravity and capillari- 
ty, the height and velocity will tends to h — H = 


20 cos 0 
: pga 
and h + 0 at t > ov, respectively. To satisfy the con- 


ditions, we can propose that the capillary rise takes the 


following form 
h(t) = H[1 — fie), (7) 


where the exponent 8 and function f(t) are to be deter- 
mined. 

Substitute the Eq.(7) and h(t) = —Hf’(t)e7® + 
H£f(t)e~® into the initial conditions h(0) = 0, h(0) = 
0, the initial condition for the function f(t) is given by 


f(0) =1, f'(0) = £. (8) 


There are many possibility on the construction of f(t) 
which can satisfy the condition in Eq.(8), for instance, 
f(t) = 1 + 6t+ at’, in which the exponent d > 2 and 
a is a constant. Since capillary rise dynamics process 
usually takes a millisecond, so the lower order of time t 
will dominate the process, based on this understanding, 
we set the function f(t) = 1+ St + atè. Hence, the 
capillary rise can be proposed as follows 


h(t) = H[1 — (1+ Bt + at*)e**]. (9) 


Thus, we have capillary velocity h(t) = H(aft? — 3at + 
B?)te-** and acceleration h(t) = —H(a6?t® + Bt — 
6aft? + bat — B?)e—**. 

The capillary rise in Eq.(9) satisfies both initial con- 
dition and infinity condition. The unknown constants 3 
can be determined by the exact series solution in Eq.(6), 


12 ocos F. Hence 
ii ? 


which gives the acceleration h(0) = 7 oa? 


by equaling h(0) = H82, we obtain the exponent 


B= 12 aøcosð — |12g 
- VTH pg V Ta’ 


To find the parameter a, we can expand the capillary 


(10) 


rise in Eq.(9) in Taylor series, and get its coefficient of 
the fifth order team of t as -+ H6?(6° + 15a), and let 
it equal to the coefficient of të in Eq.(6), namely 


864 uo? cos? 6 


1 2/93 
= = 11 
which gives 
7 0 1 12 i 
Ho cos g 
= . 12 
i paž 15 (y Ta ) (2 


It is clear to see that @ is crucial parameter which links 
to all constants the capillary rise dynamics. It is worth 


by adjusting the radius. 


to mention an interesting thing that a can be set to zero 
Finally, we obtain the capillary rise h(t) 
12 
fı- 1+ ( z) t+ at 
Ta 


1 a 
and the capillary velocity h(t) 


12 12 _(,/Bs 
(eey = 3at 4 z) te ( i e 
a a 


(14) 


20 cos 0 
h(t) _ O COS 
pga 


dh(t) _ 20 cos@ 
pga 


dt 


as well as the capillary acceleration h(t). 

The above solutions in Eqs.(13) and (14) are valid for 
the whole domain of time t € [0,00), for my knowledge, 
which have never been proposed before. 

The solution in Eq.(13) reveals that the capillary rise 
h(t) is mainly controlled by the H = ac cost with decay 


T2g 
rate of (VE), where the radius a is the only domi- 


nate parameter, the smaller radius the fast decay. 


NUMERICAL VERIFICATION AND 
DISCUSSIONS 


Verification of the solution in Eq.(13) 


To verify the correctness of the solution in Eq.(13), we 
can compare the result obtained numerically by Zhmud 
et al.[20]. for comparison study, the data of the diethyl 
either in glass shown in Table I. 


TABLE I: The diethyl ether in glass with arbitrary radius 
Ra 
2 


k k 
AEA : 9 gle) elz] 
2.2- 1074 1.67- 107? 26° 9.81 710 


ol 


In this case, for any radius a, we have the capillary rise 


H = 42..10~-°, our analytical capillary rise is given by 


4.3 -107° 
h(t) = a 
4.1 4.585- 1071! 4.597, 2, acy 
kE + ( oF “3/2? Je ve"}. (15) 


If we set the radius a = 0.0005, we have the capillary rise 
as follows 


h(t) = 0.00862 [1 — (1 + 183.397¢ + 1.056 - 10°¢?)e~ 183-3974] , 

(16) 
which is plotted in the figure 4. From this figure we can 
see that our solution agrees with [20] very well. 
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FIG. 4: Comparison study: a = 0.0005m and @ = 26° 


Capillary rise, velocity and acceleration for radius 
a = 0.0001 and wetting angle 0 = 0 


Applying all formulations to the diethyl either in glass 
with data shown in Table II. 


TABLE II: The diethyl ether in glass with radius a = 0 


ul 2) o [$8] a [m] 6 9 [5] e[2S] 
710 


2.2- 1074 1.67-1072 1x 107-4 0 9.81 


The capillary rise h(t) of the diethyl either is given by 


h(t) = 0.00479[1 — (14+ 129.68t — 94373t?)e~ 129-8], (17) 


the capillary velocity h(t) 
h(t) = —58687.26t(t + 0.027)(t — 0.05)e7 179-88, (18) 


and the capillary acceleration h(t) 

h(t) = 7.61 x 10°(¢ + 0.022) (t — 0.008) (t — 0.06)e~179-6**, 
(19) 

Those three equations are plotted in figure 5. 

Figure 5 shows that the capillary rise almost reaches 
H at time t > 0.06 seconds, and the capillary velocity 
reaches its peak point at t = 0.01 and then to zero at 
t > 0.06. The capillary acceleration reaches the peak 
point at t = 0.01 and decrease to zero at t = 0.008 and 
t > 0.06, respectively. 


Influence of radius change 


If keeping the wetting angle as 6 = 26° [20], and just 
adjusting the radius of the glass in Eq.(15). We have 
graphs in the Figure 6, which reveals that capillary rise 
is very sensitive to the radius. The smaller the radius 
in which the liquid can travel, the further up it goes. 
Therefore, the glass tube radius is a crucial parameter 
for capillary rise dynamics. The radius has also strong 
influence on the capillary velocity and acceleration which 
can been seen in the figure 7. The smaller the radius the 
fast the capillary goes. 


Discussions 


If we put all the pictures together, we can have a fair 
picture of the capillary rise dynamics as follows: the cap- 
illary rise onsets from zero initial height and velocity, at 
t = 0 as the first driving force, the surface tension sup- 
plies a kick-off acceleration in the amount of H?, this 
surface tension steadily drives the capillary rising at an 
almost uniform speed until reaches its peak at certain 
point, then the capillary speed is gradually reduced to 
zero due to the combining effect of wall viscous friction 
of galls and gravity. In the same time, the capillary accel- 
eration decrease and reaches its lowest point before the 
increase to its final static state. 
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FIG. 5: The capillary rise, velocity and acceleration for radius 
a = 0.0001 and wetting angle 0 = 0 


CONCLUSIONS 


Since the capillary rise generally has low Reynolds 
number, the flow in the tube can be considered as a lam- 


inar flow. To avoid the singularity problem, the added 


capillary rise 


FIG. 6: The capillary rise of four different radius 


mass must appear in the governing equation. From the 
dimensionless form of Bush’s equation, we found that the 
capillary rise dynamics is mainly controlled by the Bond 
number and the Galileo number. We have successfully 
obtained the perturbation and series solutions. Due to 
the poor rate of convergence of this solution, we pro- 
posed a simple approximate analytic solution by using 
series solution and initial conditions. We have verified 
the proposed analytic solution by some numerical exam- 
ples. 

Last but at least, it might be worth to point out that 
all previous solutions are for the non-oscillatory regime, 
however, when the liquid surface reaches h = H, the sur- 
face will oscillating for while before asymptotically stop. 
Within the oscillatory regime, the above solutions should 
be modified a little bit, which will be investigated by 
Homotopy Analysis Method [26, 34] in the forthcoming 
paper. 


ACKNOWLEDGEMENTS 


This paper is dedicated to the victims and heroes of 
the September 11 attacks. This work is part of the self- 
funded project: Science in Small World. 


* Electronic address: sunb@cput.ac.za 
[1] A. 


itatserscheinungen (Conclusions [drawn] from capillary 


Einstein, Folgerungen aus den  Capillar- 


velocity 


|—-m™— a=0.008 
[—— a=0.005 


acceleration 


FIG. 7: The capillary rise velocity and acceleration for two 
different radius 


phenomena), Annalen der Physik, 309 (3) : 
(1901). 
2] R. Lucas. Ueber das Zeitgesetz des kapillaren Aufstiegs 
von Flüssigkeiten. Kolloid-Zeitschrift, 23:15-22 (1918). 
3] E.W. Washburn. The dynamics of capillary flow. Physical 
Review, 17(3):273-283 (1921). 
4| E.K. Rideal. On the flow of liquids under capillary pres- 
sure. Philos. Mag. Ser. 6, 44:1152-1159 (1922). 
5| C.H. Bosanquet. On the flow of liquids into capillary 
tubes. Philos. Mag. Ser. 6, 45:525-531 (1923). 
6] W.E. Brittin. Liquid rise in a capillary tube, J. Applied 
Physics, 17:37-45 (1946). 
7] S. Newman. Kinetics of wetting of surfaces by polymer- 


513-523 


s - capillary flow, J. of Colloid and Interface Science, 
26(2)209-213 (1968). 

8] J. Szekely, A.W. Neumann and Y.K. Chuang. The rate of 
capillary penetration and the applicability of the wash- 


burn equation. J. Colloid and Interface Science, 35, 273- 
278 (1971). 
[9] J. van Brakel and P.M. Heertjes. Capillary rise in porous 


[10] 


[11] 


12 


13 


14 


15 


16 


[17] 


18 


19 


20 


21 


[22] 


[23] 


media. Nature 254, 585-586 (1975). 

S. Levine, P. Reed, E.J. Watson and G. Neale. A theory 
of the rate of rise of a liquid in a capillary. In M. Kerk- 
er, editor, J. Colloid and Interface Science, 3: 403-419. 
(Academic Press, New York, 1976). 

S. Levine, J. Lowndes, E.J. Watson and G. Neale. A theo- 
ry of capillary rise of a liquid in a vertical cylindrical tube 
and in a parallel-plate channel. J. Colloid and Interface 
Science, 73(1):136-151 (1980). 

A. Marmur. Penetration and displacement in capillary 
systems of limited size. Adv. Colloid Interface Sci., 39:13- 
33 (1992). 

M. Dreyer, A. Delgado and H.J. Rath. Fluid motion in 
capillary vanes under reduced gravity, Microgravity Sci- 
ence and Technology, 4:203-210 (1993). 

N. Ichikawa and Y. Satoda. Interface dynamics of capil- 
lary flow in a tube under negligible gravity condition. J. 
Colloid and Interface Science, 162:350-355 (1994). 

T. Delker, D.B. Pengra and P.-Z. Wong. Interface pin- 
ning and the dynamics of capillary rise in porous media. 
Physical Review Letters, 76(16):2902-2905 (1996). 

A. Marmur and R.D. Cohen. Characterization of porous 
media by the kinetics of liquid penetration: The verti- 
cal capillaries model. J. Colloid and Interface Science, 
189(2):299-304, (1997). 

R. D. Deegan, Olgica Bakajin, T. F. Dupont, G. Hu- 
ber, S. R. Nagel and T. A. Witten. Capillary flowasthe 
cause of ring stains fromdried liquid drops, Nature 389, 
827 - 829 (1997). 

D. Quere. Inertial capillarity. Europhys. Lett., 39(5):533- 
538 (1997). 

J. Fisher. Physics takes the biscuit, Nature, 397,469 
(1999). 

B.V. Zhmud, F. Tiberg, K. Hallstensson. Dynamic- 
s of capillary rise. J. Colloid and Interface Science, 
228(2):263-269 (2000). 

A. Siebold, M. Nardin, J. Schultz, A. Walliser and M. 
Oppliger. Effect of dynamic contact angle on capillary 
rise phenomena. Colloids and Surfaces A, 161(1):81-87 
(2000). 

M. Lago and M. Araujo. Capillary rise in porous media. 
Physica A, 289:1-17 (2001). 

J. Bico and D. Qur, Rise of liquids and bubbles in angu- 
lar capillary tubes. J. Colloid Interface Sci. 247:162-166 


24 


25 


26 


27 


28 


29 


30 


31 


32 


33 


34 


35 


36 


37 


38 


(2002). 

A. Hamraoui and T. Nylander. Analytical approach for 
the Lucas- Washburn equation. J. Colloid and Interface 
Science, 250:415-421 (2002). 

P. G. de Gennes, F. Brochard-Wyart and D. Quéré. Cap- 
illarity and Wetting Phenomena: Drops, Bubbles, Pearls, 
Waves. (Springer, New York, 2003). 

S. J. Liao, Beyond Perturbation: Introduction to the Ho- 
motopy Analysis Method, Chapman & Hall/CRC, Boca 
Raton, 2003. 

C. J. Lv, C. Chen, Y. C. Chuang, F. G. Tseng, Y. J. Yin, 
F. Grey, and Q. S. Zheng. Substrate curvature gradient 
drives rapid droplet motion, Phys. Rev. Lett. 113, 026101 
(2014). 

Y. Xiao, F. Yang and R. Pitchumani, A generalized anal- 
ysis of capillary flows in channels, J. Colloid and Interface 
Science, Vol. 298, pp: 880-888 (2006). 

J.W.M. Bush and D.L. Hu. Walking on Water: Bioloco- 
motion at the Interface, Annu. Rev. Fluid Mech., 38:339- 
36 (2006). 

C. Hall and W. Hoff. Rising damp: capillary rise dynam- 
ics in walls. Proc. R. Soc. A 463(2007). 

N. Fries and M. Dreyer. An analytic solution of capillary 
rise restrained by gravity, 320, 1:259-263 (2008). 

F.J. Higuera, A. Medina and A. Linan. Capillary rise of a 
liquid between two vertical plates making a small angle. 
Phys. Fluids 20, 102102, 1871 - 1884(2008) 

A. Ponomarenko, D. Quéré and C. Clanet. A univer- 
sal law for capillary rise in corners, J. Fluid Mech., vol. 
666:146-154 (2011). 

S. J. Liao, Homotopy Analysis Method in Nonlinear Dif- 
ferential Equations, Springer-Verlag, New York, 2011. 
Y.P. Zhao. Physical Mechanics of Surface and Interface 
(in Chinese) (Science Press, Beijing, 2012). 

W.M. Bush. Interfacial Phenomena, MIT Open Course 
(2013). 

O. Shardt, P. R. Waghmare, J.J. Derksena and S. K. 
Mitrab. Inertial rise in short capillary tubes, RSC Adv., 
4:14781-14785 (2014). 

X.X. Dou, S.P. Li, and J.L. Liu. Zero curvature-surface 
driven small objects, Appl. Phys. Lett. 111, 081602 
(2017). 


